Cavitation pressure in liquid helium 
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Recent experiments have suggested that, at low enough temperature, the homogeneous nucleation 
of bubbles occurs in liquid helium near the calculated spinodal limit. This was done in pure superfluid 
helium 4 and in pure normal liquid helium 3. However, in such experiments, where the negative 
pressure is produced by focusing an acoustic wave in the bulk liquid, the local amplitude of the 
instantaneous pressure or density is not directly measurable. In this article, we present a series of 
measurements as a function of the static pressure in the experimental cell. They allowed us to obtain 
an upper bound for the cavitation pressure Pcav (at low temperature, Pcav < —2.4 bar in helium 3, 
Pcav < —8.0 bar in helium 4). From a more precise study of the acoustic transducer characteristics, 
we also obtained a lower bound (at low temperature, Pcav > —3.0 bar in helium 3, Pcav > —10.4 bar 
in helium 4). In this article we thus present quantitative evidence that cavitation occurs at low 
temperature near the calculated spinodal limit (—3.1 bar in helium 3 and —9.5 bar in helium 4). 
Further information is also obtained on the comparison between the two helium isotopes. We finally 
discuss the magnitude of nonlinear effects in the focusing of a sound wave in liquid helium, where 
the pressure dependence of the compressibility is large. 



PACS Numbers: 67.55. Cx, 67.40.Kh, 64.60. Qb, 43.25.+y 
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I. INTRODUCTION 

As has often been explained, the purity of liquid heliurn. 
offers unique opportunities to study phase transitions Jj'EI 
Furthermore, liquid helium is a simple system and its 
thermodynamic properties hfV^ been accurately mea- 
sured, so that several authorsQTl^ have calculated the ex- 
trapolation at negative pressure of its equation of state. 
All these calculations are reasonably consistent with each 
other; they predict the existence of a spinodal limit Pg 
at —3.1 bar for liquid helium 3 and_at —9.5 bar for liquid 
helium 4, near the absolute zero.Ej As one approaches 
the spinodal limit, the compressibility diverges so that 
the energy barrier for the nucleation of bubbles in the 
stressed liquid vanishes, and the liquid becomes totally 
unstable. 

In order to test this prediction, we have studied cavi- 
tation by focusing a high intensity acoustic wave in bulk 
liquid helium, far from any wall. If one wants to approach 
the spinodal limit, it is necessary to work at very low tem- 
perature, otherwise thermal fluctuations allow the system 
to pass a rather high energy barrier at higher, i.e., less 
negative, pressure. This is possible in liquid helium only, 
since no other liquid exists down to zero temperature. 
Eventually, close enough to the spinodal limit and at low 
enough temperature, a crossover has been predicted to 
exist from a thermally activated classical cavitation to 
a quantum regime where the nucleation of bubbles oc- 
curs twptiie quantum tunneling of a sizable quantity of 
liquidfiy 

In previous experiments, we obtained some evidence 
for the existence of this crossover in helium 4. Below 
about 0.6 K, Lambare et al&3 observed a stochastic and 
temperature independent behavior, while above 0.6 K 



they observed another stochastic behavior where the cav- 
itation threshold decreased with temperature. These 
observations agreed with the theoretical predictions by 
MarisoQ if the quantum cavitation occurred at —9.2 bar, 
close enough to his calculated value of the spinodal limit. 
This agreement was obtained after noticing that, given 
the amplitude of the wave, and due to adiabatic cooling 
during the negative pressure swing, the minimum instan- 
taneous temperature was lower than the static temper- 
ature by a factor of 3 (that is, 0.2 K for a static tem- 
perature of 0.6 K). Unfortunately, in this experiment, it 
was not possible to measure directly the exact value of 
the negative swing in the acoustic wave when nucleation 
occurs. Instead we measured the voltage applied to the 
transducer that generates the wave; as for the wave am- 
plitude in the focal region, it is difficult to estimate since 
one expects the focusing of the wave to be nonlinear and 
its precise calculation in a cylindrical geometry has not 
yet been done. 

Given these results and difliculties, we tried a compari- 
son with helium 3. We first confirmed that the cavitation 
threshold pressure is less negative in helium 3 than in he- 
lium 4.II3 Moreover, we did not observe a crossover to a 
temperature independent quantum regime of cavitation. 
The latter observation may be understood_bji invoking 
theoretical arguniflBts involving dissipationlij'tZI or Fermi 
liquid properties Ji'tS 

In this context, we have tried to improve our estimation 
of the cavitation pressure by carefully measuring the de- 
pendence of the cavitation voltage on the static pressure 
Pstat in the cell. As explained in this article, this study 
allows us to present an upper bound for the cavitation 
pressure in both helium 3 and helium 4 (Pcav < —2.4 bar 
in helium 3, Pcav < —8.0 bar in helium 4). From a more 
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FIG. 1. Sketch of the experimental setup: a high amplitude 
pressure swing is generated at the focus of the hemispherical 
transducer; a laser beam passes through the acoustic focal re- 
gion, and the light scattered by cavitation is detected with a 
photomultiplier tube (PMT). 

precise study of the acoustic transducer characteristics, 
we also obtained a lower bound (Pcav > —3.0 bar in he- 
lium 3, Pcav > — 10.4bar in helium 4). Finally, we also 
encountered interesting questions about the magnitude of 
nonlinear effects in the focusing of high intensity acous- 
tic waves in a medium such as liquid helium where the 
compressibility strongly depends on pressure. 
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FIG. 2. Cavitation probability versus driving voltage in he- 
lium 3 at T = 566 mK and P = 4 mbar. Each circle is a prob- 
ability measurement over 1600 bursts at a given voltage. The 
solid line is a fit with Eq. (jl|) which gives 14 = 22.403 V and 
5 = 355. 

tube position and the size of a diaphragm in front of it, 
we can choose to detect cavitation either from the scat- 
tered light or from the light missing in the transmitted 
beam. 



II. EXPERIMENTAL PROCEDURE 

A. Experimental setup 

As shown in Fig. ^ our experimental, satup-jis similar 
to the one used in previous experiments.EJIl3'E3 We have 
reduced the inner volume of the cell to 4.5 cm"^ in order to 
lower the cost of the helium 3 experiment. This cell is an- 
chored to the mixing chamber of a dilution refrigerator 
by copper columns with gold plated contacts. Carbon 
resistor thermometers and a heater are thermally con- 
nected to the outside of the copper cell walls. A burst of 
1 MHz ultrasound is emitted and focused in the liquid by 
a hemispherical piezoelectric transducer (the "ceramic" ) . 
The cell is connected to a buffer volume at room temper- 
ature, so that we easily monitored the static pressure in 
the cell. When working with helium 3, this was done 
with a capacitive sensor (Keller type PAA-41) with an 
accuracy of ±0.5 mbar. For helium 4, we used conven- 
tional pressure gauges; the uncertainty was ±0.7 mbar 
for low pressure measurements and ±12.5 mbar for static 
pressures above Ibar. Two sets of four windows allow 
us to shine a laser beam from the outside of the cryostat 
through the acoustic focal region, and to detect the light 
that is scattered by local changes in density. The detec- 
tor is a photomultiplier tube with a response time shorter 
than 0.5 /iS. The acoustic wave scatters the laser light at 
small angle; when cavitation occurs, gas bubbles scatter 
light at a larger angle. By adjusting the photomultiplier 



B. Voltage measurements 

One side of the ceramic is grounded in the cell and 
the other is connected to the output of a rf amplifier 
at room temperature. In order to improve the accuracy 
of the probability measurements, we needed to increase 
the number of events to count. Furthermore, to study 
helium 3 down to 40 mK, we had to reduce the dissi- 
pation by lowering the repetition rate of the acoustic 
pulses. Indeed, inside the ceramic itself, there is a me- 
chanical dissipation of order 2.6 /j,W for six-cycle pulses 
repeated at IHz, a non-negligible amount. As a conse- 
quence, we had to improve the stability of the excitation. 
Since this piezoelectric transducer has a low impedance 
at resonance (12 to 16 fl, see below), we had to mini- 
mize the possible drift of the impedance of its connect- 
ing cable. The previous cable included a 1.5 m section 
from 4K to room temperature, that was a commercial 
cable made of stainless steel. We replaced this section 
by a homemade coaxial cable consisting of a 0.24 mm di- 
ameter copper wire separated from stainless steel tubing 
(1.2mmi.d., 1.5mmo.d.) by a Teflon tube (0.30mmi.d., 
0.76mm o.d.). We kept the superconducting coaxial ca- 
ble connecting the transducer to the 4K region. The 
impedance of the whole line is now negligible compared 
to the transducer impedance, so that the excitation does 
not depend on the level of helium in the 4K bath. We 
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finally improved the stability of the excitation voltage it- 
self by using a new generator (Hewlett Packard model 
33 120 A). Our homemade rf amplifier uses an Apex PA09 
circuit |-and is located in a thermally regulated box as 
before JlJ 

During the experiment, the ceramic is driven by a burst 
with a few cycles of a 1 MHz sine wave. This burst is mon- 
itored with a digital oscilloscope (Tektronik model TDS 
420A, 200 MHz, 100 MS s"!). The waveform is trans- 
ferred to a computer using a conventional IEEE inter- 
face and then fittecUtp a sine wave using the Levenberg- 
Marquardt methodcJ with four adjustable parameters 
(amplitude, frequency, phase, and offset). The first cycle 
is always slightly distorted and therefore ignored in the 
fit. We have checked the distribution of amplitudes over 
10^ fits under the same experimental conditions. The 
repetition rate being 1 Hz in this case, the total acqui- 
sition time was about 28 h. Over such a duration, we 
achieved a stability of 1.8 x 10^"^. 



C. Statistics of cavitation 

At a given temperature and pressure in the cell, 
we measured the cavitation probability S over series 
of bursts_.ai. several driving voltages. As previously 
reportedoli3 and illustrated in Fig. ||, cavitation is a 
stochastic process, and S(V) is well described by the 
"asymmetric S-curve formula" 



E(y) = 1 - cxp 



-ein2exp( ^-1 



(1) 



We previously explained^ that this behavior is ex- 
pected if cavitation is thermally activated. It is derived 
from a linear expansion of the activation energy E around 
the cavitation threshold voltage Vc where the probability 
is 1/2. The quantity ^ is the inverse width of the dis- 
tribution of events, given by ^ = {Vc/kBT){dE/dV). A 
similar formula is also expected in the quantum regime: 
^ = {Vc/K){dB / dV) where B is an action. In our whole 
study, the cavitation threshold Vc is determined from a 
fit of S with Eq. (|l|); the accuracy is ±0.1 %. 



D. Piezoelectric transducer characteristics 

In the following discussion, we will need the various 
parameters that describe the transducer behavior. It is 
made of lead zirconium titanate (Quartz & Silice type 
P7-62); its thickness is 2 mm, its inner radius i?trans is 
8 mm, and its mass M is 7.5 g. It is used at the reso- 
nance of its first thickness mode (/o = 1022 kHz). Near 
its resonance, wc found that its inverse impedance or ad- 
mittance y = 1/Z is described by a Lorentzian curve 
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where ujq = 27r/o. After measuring i?, we fitted the above 
equation to obtain the quality factor Q in helium 3 or 4 
at low temperature. We obtained the preliminary values 
g = 175 ± 25 for helium 3 and Q = 135 ± 10 for he- 
lium 4. Unfortunately, the thickness mode resonance is 
close to another resonance which must be a high order 
fiexion mode, so that we do not take this determination 
of Q as reliable. 

The quality factor is of particular interest for the short 
sinusoidal pulses that we are using. Indeed, it character- 
izes the transient time that is necessary to build up the 
amplitude of the oscillation of the transducer. At a fre- 
quency /o, this transient time is tq = Q/(27r/o). For the 
sake of simplicity, we will take the period as the time 
unit; in other words, a time t corresponds to x cycles 
through the relation x — fot. Let Vq be the amplitude of 
the sinusoidal voltage and n the total number of cycles 
in the electrical burst; hence the burst duration is n in 
our reduced units. During the pulse (0 < a; < n), the 
intantaneous voltage is V{x) = Vb sin(27ra:). 

In a steady regime, i.e., in the long time limit, the 
ceramic surface would oscillate with an amplitude 
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This equation is derivedcil by writing 
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where -Egtorod is the energy stored during one period 
and Pdiss = Vb^/(2_R) the average dissipated power. 
Let / be the transducer thickness and C(^7^) the ampli- 
tude of the displacement at a distance z from the in- 
ner surface. For the lowest thickness mode, which is 
a standing wave with one node at z — 1/2, we have 
C(z,t) = ^oCOs(7rz/Z) sin(wot). The density of acous- 
tic energy is the sum of a kinetic term and a pressure 
term that have the same amplitude. By integrating this 
energy density over the transducer volume and averag- 
ing over one period, we find the total acoustic energy 
^-acoust = Mloq'^^o'^ /A. For a piezoelectric transducer, 
there is an additional term in the stored energy: 
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where C is the transducer capacitance and kt the elec- 
tromechanical coupling constant. We measured C = 1 nF 
at low temperature. As for fct, the constructor gives 0.47 
at room temperature, aad an experimental check at low 
temperature gave 0.27.Eil With these values and the mea- 
sured resistance at resonance i? = 12 17, we find that the 
ratio between the second term in Eq. (||) and -Estorod cal- 
culated from Eq. (^) is of the order of 5 x 10~^. There- 
fore, we can neglect the second term of Eq. (|^) and write 
E^acoust = E^storod to derive Eq. (|). 

In the transient regime, x cycles after the excitation 
has started, the actual surface oscillation is 
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FIG. 3. Cavitation threshold voltage versus number of cy- 
cles in helium 3 (circles, T = 423 mK and P = 4mbar) and 
helium 4 (squares, T = 284 mK and P = 40mbar). The solid 
lines are fits with Eq. (^ which gives Q = f 19 in helium 3 
and Q = 100 in helium 4. 



C(x) = Co sin(27ra;) 
for < a; < n, and 
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C(x) = Co sin(27ra;) 



1 — exp 



Q 



exp 



-27r- 



Q 



(7) 



for X > n. This leads to a better method for the measure- 
ment of the quality factor. We have used the dependence 
of the cavitation threshold voltage on the total number 
of cycles n in the electrical pulse. 

Indeed, cavitation occurs for a well defined value of 
the surface oscillation corresponding to a well defined 
pressure at the acoustic focus. It is the most negative 
pressure in the acoustic wave. Given the phase of the 
excitation voltage at time t = 0, cavitation occurs for 
Xc = n + J, and we compared our measurements of the 
cavitation voltage Vc (n) to the equation 
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As shown on Fig. |^, we obtained excellent fits with 
Eq. (^. In helium 3, we found (3 = 119±4. In hehum 4 
we found Q — 100 ± 2. In the same run, we also veri- 
fied that the resonance frequency (1022 kHz) minimizes 
Vc- We will use these values of Q in the analysis below. 
We limited the fit to n < 21 because for Q = 100 and 
n > 23 cavitation occurs during the previous swing, at 
time Xc = n — J. 

In principle, one expects slightly different values for 
Q in helium 3 and in helium 4. The coupling of the 



ceramic to liquid helium is small, because the acoustic 
impedance pc in liquid helium (of order 10^m~^s~^) is 
small compared to the impedance pc Cc of the ceramic 
(3 X 10^ m-2 s-i). The ratio of the emitted power Pa to 
the dissipated power Pq in the ceramic is given by 



P. ^ 4 pc 



(9) 



As a consequence, the emitted acoustic power should be 
6% for helium 3 and 15% for helium 4, in qualitative 
agreement with our measurement. As for the electrical 
impedance at resonance, we found 12.2 Q in helium 3 and 
15.6 Q in helium 4. 

Having chosen the values of the frequency (1022 kHz) 
and the pulse width (three to six cycles), we studied the 
temperature and pressure dependence of the cavitation 
threshold. 



E. Thermal relaxation 

To ensure reproducibility of the measurements, we had 
to check the time needed by the system to reach a new 
steady state after each adjustment of the temperature 
regulation. In order to do this, we used the following 
procedure. Starting from a low temperature (65 mK, for 
instance) , we set the driving voltage to a value such that 
the probability was around 20%; then we warmed up 
to a new temperature (10 mK higher, for instance): the 
reading of our thermometers was almost instantaneous, 
whereas the probability increased slowly. This is because 
small temperature gradients between the acoustic focal 
region and the thermometer take some time to vanish 
completely (see Sec. HF ). Then we adjusted the driving 
voltage to keep the probability around 50%, where the 
system is most sensitive to temperature drifts. We waited 
until the probability measurement fluctuated around a 
constant value. After several checks around different 
temperatures, we found the relaxation time to be less 
than 90 min at low temperature in helium 3; we thus de- 
cided to wait for 90 min after each temperature change in 
the lowest temperature range. The relaxation time was 
found to be less than 30 min for helium 3 above 200 mK 
and less than 5 min for helium 4 at all temperatures, so 
that it was sometimes possible to work faster. In ad- 
dition, we always waited several hours after each liquid 
helium or nitrogen transfer, because longer drifts were 
observed, particularly in the optical detection of the bub- 
bles. The latter problems are associated with the thermal 
contraction of the refrigerator. 



F. Dissipation at low temperatures 

Working in helium 3 at low temperatures raises prob- 
lems of temperature inhomogeneities. Unlike superfluid 
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FIG. 4. Cavitation threshold voltage versus temperature: 
(a) in liquid helium 3 (P = 12mbar, 3 cycles); (b) in liquid 
helium 4 (P = 37mbar, 6 cycles). Empty circles represent 
raw data, and full circles data corrected from the sound at- 
tenuation. 

helium 4, helium 3 is a poor thermal conductor. There- 
fore, the dissipation in the cell leads to temperature gra- 
dients between the cavitation area and the cell walls 
where the temperature is measured. There are two main 
sources of dissipation: mechanical friction in the ceramic 
and absorption of the laser light in the cell. The former 
was minimized by decreasing the pulse duration and the 
repetition rate; the latter by attenuating the laser inten- 
sity and chopping the beam. We could estimate both ef- 
fects by monitoring the heating power of the temperature 
regulation under different conditions. For the mechanical 
dissipation, we found 2.6 /iW for six-cycle pulses repeated 
at 1 Hz. We then chose to work with the shortest pulses 
allowed by the power of the amplifier and the lowest af- 
fordable repetition rate, namely, three cycles and 0.05 Hz. 
We checked that the cavitation voltage was independent 
of the repetition rate below 0.1 Hz. 

As for the light absorption we found it to be 0.7 fiW in 
a cw operation mode, after having attenuated the laser 



beam. For a further reduction, we built an optical chop- 
per that was synchronized with the acoustic pulse. With 
an ordinary electrical relay, we easily achieved an opening 
time of 20 ms, so that the radiation power was divided 
by a factor of 1000 for a repetition rate of 0.05 Hz and 
became completely negligible. 

In helium 4 these problems did not exist and we worked 
with repetition rates ranging from 0.5 to 2 Hz and a pulse 
width of six cycles. The lowest rates were used to reach 
the lowest temperatures; we kept the optical chopper for 
the same reason. We had to use longer pulses because 
cavitation requires a larger amplitude than in helium 3, 
and our rf amplifier has a limited output amplitude. 



G. Sound attenuation 



In all our experiments, we had to consider the atten- 
uation of sound. In helium 3, this attenuation becomes 
important when the liquid enters the Fermi liquid region; 
this is below 100 mK, where the viscosity varies as 1/T^ 
because of the-temperature variation of the quasiparticle 
collision time.c3 We used_the values of the viscosity mea- 
sured by Bertinat et a/.Ej and the best fits they give in 
two temperature regions (Eqs. (7) and (13) of Ref. |23| ). 

In helium 4, the temperature variation of the sound 
attenuation shows a peak at the temperature where 
the phonon-roton collision time eauals the sound period 
(750 mK for a 1 MHz sound wave) rH We used the absorp- 
tion at 1 MHz extrapolated by Maria^i from the measure- 
ments of Waters et alE5 at higher frequency. 

In both cases, we assumed that, outside the focal re- 
gion whose typical radius is one acoustic wavelength, the 
sound wave has a small enough amplitude to be treated 
as linear. The measured cavitation threshold voltage be- 
ing related to the emitted amplitude, we assumed that 
the amplitude at the acoustic focus is reduced by the at- 
tenuation factor exp(— aiJtran), where i?tran — 8 mm is 
the transducer radius. We finally need to remark that in 
helium 4 we had checked this attenuation in our previ- 
ous experiment, but in helium 3 a similar check appeared 
impossible to do. In helium 4, the value of the attenu- 
ation had been checkedtil by measuring the amplitude 
of the light scattered by the acoustic wave. The tem- 
perature variation of this "acoustic signal" was found in 
good agreement with the known sound attenuation up 
to 0.9 K; however, above 0.9 K, the sound amplitude was 
found smaller than predicted by the sound attenuation 
only, as if additional mechanisms had to be considered 
also (diffraction, temperature drift of the transducer reso- 
nance frequency, etc.). In helium 3, a similar check would 
have required too much light amplitude and an average 
over too many bursts to be made in a reliable way in 
the interesting temperature region, which is in the low 
temperature limit. 
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FIG. 5. Corrected cavitation threshold voltage versus tem- 
perature at different static pressures: (a) in liquid helium 3 
(3 cycles); (b) in liquid helium 4 (6 cycles). 

H. Data corrections 

Let us start with helium 4. Figure ^(b) presents a set 
of results for the cavitation threshold voltage as a func- 
tion of the static temperature T in the cell. The results 
obtained confirm what was already published by Lam- 
bare et aZ.Ej There are two sets of data (raw data and 
corrected data). Fitting cavitation probability measure- 
ments with Eq. (|l|) gives the raw value of the threshold 
voltage Vc\ it is the excitation amplitude at which the 
probability is 0.5. The correction accounts for the at- 
tenuation of the acoustic wave over its flight distance, 
i.e., the transducer radius i?tran- The corrected data in- 
dicate the existence of a crossover around T = 0.6 K. 
Below 0.6 K, we attribute the temperature independent 
regime to quantum cavitation occurring by quantum tun- 
neling. Above 0.6 K, the cavitation threshold decreases 
with increasing T, as expected for a thermally activated, 
classical regime. Our goal is to show that the cavitation 
threshold in the low temperature limit corresponds 



FIG. 6. Static pressure as a function of the parameter pVc 
(see text): (a) in liquid helium 3 (3 cycles); (b) in liquid he- 
lium 4 (6 cycles). The pressures obtained by linear extrapo- 
lation are given in Table |. 

to a pressure close to the spinodal limit Pg ~ —9.5 bar. 
As explained in Ref. |l^, agreement with the quantum 
cavitation theory of Maris requires that this quantum 
plateau corresponds to —9.2 bar. We have tried to check 
this value more precisely in this work. 

As for helium 3, our results are presented in Fig. ^a). 
We have extended to lowa^temperatures the measure- 
ments we reported before.liSi A striking difference from 
helium 4 is the sharp increase of the threshold in the 
very low temperature limit (below 55 mK). It occurs in 
a temperature domain where the applied correction di- 
verges. Indeed, the attenuation of sound is proportional 
to in the Fermi liquid region. One possible origin 
of the experimental divergence could be that we do not 
apply the right correction. However, we do not believe 
so. If this effect is not an experimental artifact, it is 
possible to propose an interpretation for it by consider- 
ing the particwlar properties of a Fermi liquid near its 
spinodal limitJla In this article, we restrict ourselves to 
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the temperature domain above 55 mK with the follow- 
ing question: does the cavitation pressure approach the 
calculated spinodal pressure Pg = —3.15 bar there? 



III. MEASUREMENTS UNDER PRESSURE 

We now turn to the dependence on the static pressure 
in the cell. If there was no attenuation of sound in liq- 
uid helium, if the coupling of the ceramic to the helium 
was independent of pressure, and if the sound amplitude 
at the focus was simply proportional to the driving volt- 
age, then the estimation of the cavitation pressure would 
be very simple. By measuring the cavitation threshold 
voltage Vc as a function of the static pressure, we would 
observe a linear variation and then extrapolate down in 
pressure. The negative pressure at which I4 = would 
be the pressure at which cavitation occurs in our exper- 
iment. Although things are not that simple, mainly be- 
cause, in helium, the sound velocity strongly depends 
on pressure so that the focusing of sound is nonlinear, 
we have used the pressure dependence of the cavitation 
threshold as explained below. 



A. Results 



are small. From the measurements of Dransfeld et al. 
we estimated the relative variation of the viscosity as 
— 2.4%bar~"'^ and that of the attenuation peak tempera- 
ture as -1-1.8 % bar~^ . Accordingly, we used a single value 
for the viscosity rj in our whole study. 

Using Ref. and references therein, we checked the 
temperature dependence of the density and of the sound 
velocity in our temperature ranges. The temperature 
variation of the density is smaller than 0.4 % in helium 3 
and smaller than 0.03% in helium 4. The temperature 
variation of c in helium 4 is less than 0.2 %. As a conse- 
quence we neglected these temperature variations in our 
analysis (for more details, see Sec. HID). 

The only significant variations are the pressure varia- 
tion of the density and that of the sound velocity. We 
took them into account in deriving the absorption coef- 
ficient: 



a 



3 pc^ 



(10) 



For the calculation of a, we used the equation of state by 
Maris. For the present purpose, it would make no signif- 
icant difference to use the one given in the Appendix. 

We can now explain how we used the dependence of 
the cavitation voltage on the static pressure to measure 
the cavitation pressure. 



We repeated our measurements at different pressures 
using an identical set of temperatures for each pressure. 
We were limited to a few bars only because, as the 
pressure increases, bubbles become more and more dif- 
ficult to detect- Indeed, as expected from the Rayleigiu 
Plesset theoryBa and previously measured in helium 4,EII 
the maximum radius of bubbles varies as eo^^'^Pstat""'^^^, 
where eg is the energy acquired by the bubble from the 
acoustic wave. In helium 3, there is less acoustic energy 
available to make the bubble grow after its nucleation. 
Above 0.5 bar in helium 3 and 2 bar in helium 4, the 
bubbles are too small to be detected at the cavitation 
threshold. Our results are shown in Fig. |[ 

A correction has been applied to these data, which ac- 
counts for the sound attenuation. In order to do this, we 
first noticed that, in both helium 3 and helium 4, the vis- 
cosity does not significantly vary with pressure from to 
3 bar. For example, in helium 3, the viscosity 77 varies as 
T^^ and the low temperature liroit of 77 varies by less 
than 6 % in this pressure rangeEj Furthermore, in this 
pressure study, there is only one measurement (namely, 
at 57 mK) for which attenuation is not negligible; we dis- 
carded measurements below 55 mK in helium 3 because 
of their very large dependence on temperature which in- 
duces some scatter in the results and makes the pressure 
extrapolation imprecise. 

In helium 4, the pressure dependence of the attenua- 
tion is not well known. Several groups have shown that 
the reduced viscosity decreases and the peak temperature 
increases with increasing pressure. However, these effect 



B. An upper bound for the cavitation pressure 

If the focusing of the sound wave was linear, the pres- 
sure swing at the focus would be 

AP = pC^2^tra„C, (11) 

where ^ is the displacement of the transducer wall. This 
displacement itself is proportional to the voltage applied 
to the ceramic since the ceramic oscillations are small 
enough to be in the linear regime. As a consequence, 
we have plotted the static pressure as a function of the 
product pVc] again, to calculate p for the different pres- 
sures, we used the equation of state of Maris and we 
neglected the temperature variation. The result is shown 
in Fig. ^(a) for helium 3 and Fig. |^(b) for helium 4. We 
then extrapolated our measurements linearly to pVc = 0. 
We claim that this extrapolation gives an upper bound 
for the cavitation pressure. 

Of course, if we knew the nonlinear relation between 
the pressure oscillation at the focus and the driving volt- 
age V, we could directly obtain the value of the cavitation 
pressure; but we do not know this nonlinear relation: we 
only guessed it as shown by the solid line in Fig. 0. One 
thing we know is the qualitative effect of nonlinearities. 
Suppose that one could work with a negative static pres- 
sure in the cell, close to the cavitation pressure. A small 
amplitude sound wave would then be sufficient to pro- 
duce cavitation, and the focusing of this small amplitude 
wave would be in its linear regime, as described by 
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FIG. 7. Illustration of the method used to obtain bounds 
for the cavitation pressure (see text). For the sake of clarity, 
this figure is not drawn on scale. 



effects are larger and larger. As a consequence, we mea- 
sure a slope dPstat/d{pVc) that is smaller than if we had 
no nonlinear effects; this corresponds to the negative cur- 
vature of the solid line in Fig. ^. A linear extrapolation 
thus leads to a pressure less negative than the actual 
pressure at the focus; this means that an upper bound 
for Pcav is obtained. This effect of nonlinearities has been 
confirmed by the preliminary resiilts of the recent numer- 
ical calculations by^Appert et oZ.E3 and by a very recent 
experimental studyEil As shown in Fig. ^ and in Table |, 
we repeated this extrapolation for series of data at differ- 
ent temperatures; at low temperature, we found —2.4 bar 
for helium 3 and —8 bar for helium 4. These two bounds 
are only slightly larger than the calculated spinodal lim- 
its of these two liquids (—3.15 and —9.65 bar), in very 
good agreement with our expectations. 



C. A lower bound for the cavitation pressure 



Eq. ([Til). Since we cannot start from a negative static 
pressure, we are forced to use a large amplitude sound 
wave. To make a given pressure swing we need to drive 
the transducer with a larger voltage than if the regime 
were linear. This is because, as one reaches more and 
more negative pressures at the focus, the liquid there is 
more and more compressible, so that the excitation is 
less and less efficient in building up a negative pressure 
swing. If we considered the effects of nonlinearities on the 
magnitude of the positive pressure swing, it would be the 
opposite. We thus expect the large amplitude sound os- 
cillation at the focus to be asymmetric, that is, to have 
a smaller negative swing amplitude and a larger positive 
swing amplitude than a sine wave. As we increase the 
static pressure, we need an even larger sound amplitude 
to reach the cavitation threshold, so that the nonlinear 



The same reasoning also leads to a lower bound for 
the cavitation pressure. Indeed, we can try to estimate 
the displacement C, and use Eq. ( pT| ) to obtain the mag- 
nitude of the pressure swing at the focus; the calcu- 
lated value will be an overestimate of the actual magni- 
tude. According to Eq. (||), when starting from a static 
pressure Pstat and a temperature T, the pressure swing 
AP(Pstat,!r) required to produce cavitation is propor- 
tional to p(P,tat,T)K°°(^'stat,T), where K°°(Pstat,T) is 
the cavitation voltage extrapolated to infinite pulse du- 
ration using Eq. (||). More precisely, we write 



AP(P,tat,T) 



2 Ptr 



/0(ntat,T)F- (Pstat, T) 

In the linear approximation, we have 



MR 



(12) 



TABLE I. Upper and lower bounds for the cavitation pressure in helium 3 and helium 4 at 
several temperatures. These values are used in Fig. H 





Helium 3 






Helium 4 




Temperature 


p 

max 


mm 


Temperature 


p 

max 


p . 

mm 


(mK) 


(bar) 


(bar) 


(mK) 


(bar) 


(bar) 


56.8 


-2.39 


-3.00 


49.8 


-8.06 


-10.35 


234 


-2.36 


-2.86 


129 


-8.06 


-10.37 


418 


-2.25 


-2.78 


216 


-7.98 


-10.44 


555 


-2.19 


-2.71 


292 


-7.98 


-10.44 


768 


-2.13 


-2.60 


414 


-8.08 


-10.40 


1085 


-1.96 


-2.37 


525 


-8.21 


-10.38 








614 


-8.28 


-10.41 








652 


-8.22 


-10.38 








702 


-8.14 


-10.21 








749 


-7.89 


-10.02 








804 


-7.77 


-9.95 








854 


-7.73 


-9.93 








901 


-7.69 


-9.81 
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AP(P,tat,T) = AP(0,r) + P,, 



(13) 
with 



so that Pstat is a hnear function of the product pV° 
a slope given by Eq. jl^). 

By inserting our measured values of R and Q 
in Eq. (p^, we estimate the slope of the linear 
regime to be 1035 Pa kg^^ V^^ m'^ in helium 3 and 
899 Pa kg^^ V^^ m'^ in helium 4. Now drawing a straight 
line with this calculated slope through the low static pres- 
sure data points, we find the lower bound for Pcav as the 
intersection with the vertical axis 

= 0); Fig. 0il- 

lustrates this construction. We repeated this procedure 
for series of data at different temperatures; at low tem- 
perature, we found —3.0 bar in helium 3 and —10.4 bar 
in helium 4. 

The bounds obtained with this second method are 



more negative than the bounds obtained in Sec. IIIB; 
this is equivalent to the fact that the experimental slopes 
of Sec. [II B, after extrapolation to infinite pulse dura- 
tion, are smaller (845 ± 16 Pa kg^^ m'^ in helium 3 
and 660 ± 15 Pa kg^^ V^^ m'^ in helium 4) than the ones 
we calculated here. Given the experimental difficulties 
that we mentioned above, this is rather satisfactory and 
it supports our whole analysis. 



D. Discussion 

As now shown in Fig. |^, we have obtained experimen- 
tal bounds for the cavitation pressure in helium 3 (a) 
and in helium 4 (b). We have estimated the error bars 
on the data points of Table || and Fig. ||. Let us start with 
Pmax- The uncertainty is about ±0.05 bar for helium 3 
and ±0.2 bar for helium 4. For helium 3, we have tried 
to include the temperature dependence of the density by 
using KoUar and VoUhardt's analysis and program;E3 the 
effect is to lower all the points by less than 30 mbar. As 
for Pmin, the main uncertainty comes from the deter- 
mination of the quantities Q and R, which are used to 
calculate the slope from Eq. ( pi] ) . This may lead to a sys- 
tematic error of about ±5%, i.e., ±0.15 bar in helium 3 
and ±0.5 bar in helium 4. In the determination of Pmin, 
there is also a small uncertainty of ±0.2 % from the mea- 
surement of pVc- 

Within about 10%, our results agree with the calcu- 
lated spinodal limits at low temperature. Of course, we 
confirm that helium 3 is about three times more fragile 
than helium 4, in the sense that cavitation occurs (the 
liquid breaks) at a pressure that is about three times less 
negative in helium 3 than in helium 4. In fact, what sur- 
prises us on this figure is the smallness of the difference 
between the upper bounds and the lower bounds. This 
difference is related to the magnitude of the nonlinear 
effects, and we find these nonlinear effects rather small. 
The smallness of the nonlinearities can also be seen in the 
fact that our measurements at different temperatures fall 
on parallel lines: the slope of Pstat (pK) is nearly constant 
although Vc varies. This is a somewhat surprising 
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FIG. 8. Bounds obtained for the cavitation pressure in he- 
lium 3 (a) and helium 4 (b). Upper bounds are given by the 
circles and lower bounds by the squares. The arrows indicate 
the theoretical values of the cavitation press ure at low temper- 
ature. Error bars are discussed in Sec. |IIID| ; their magnitudes 
for the lower bounds are smaller than the marker size. 

behavior that would need to be confirmed by a direct 
measurement of the sound amplitude at the focus. How- 
ever, with such an open geometry where there is nothing 
in the acoustic focal region, the optical measurement of 
the instantaneous sound amplitude seems very difficult 
to perforfHi, as shown by the previous attempt by Nis- 
sen et alSn It would thus be very interesting to calculate 
these effects in our case of a hemispherical transducer. 
The only calculations performed up to now are done in 
a fully spherical geometry because it is one dimensional 
(everything depends on the distance r to the center only). 
According to the existing calculations if-ithis spherical ge- 
ometry, the nonlinear effects are large.Ej A calculation in 
a hemispherical geometry looks much more difficult be- 
cause it is two dimensional. Its results might be different 
because the local condition at the center is not the same: 
there is a velocity node at the center in the spherical 
geometry which is closed, and not in the hemispherical 
geometry which is open. 
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IV. CONCLUSION 

By studying the pressure dependence of cavitation in 
liquid helium, we have obtained bounds for the cavitation 
pressure: at low temperature, —3.0 < Pcav < —2.4 bar 
in helium 3 and —10.4 < Pcav < —8.0 bar in helium 4. 
These negative pressures are close to the calculated spin- 
odal limits (—3.1 bar in helium 3 and —9.5 bar in he- 
lium 4), as expected for homogeneous nucleation near 
absolute zero. In order to improve this accuracy, we be- 
lieve that it is necessary to insert something in the cavita- 
tion region, for example a glass plate. However, in such a 
case, the plate may affect the cavitation conditions—This 
type of experiment is in progress in our laboratory.Ld Our 
measurements also give the temperature dependence of 
the cavitation pressure: we will discuss it in a forthcom- 
ing paper. 
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APPENDIX: 

The method used by MarisBi to obtain the value of the 
spinodal pressure Pg consists in extrapolating measure- 
ments of the sound velocity c at positive pressure with a 
law of the form 

c^ = b{P-Ps). (Al) 

Using the measurements of Abraham et aZ.js^S Maris 
found the value of Pg to be —3.097 bar in helium 3 (Ref. |^) 
and —9.5219 bar in helium 4 (Ref. ||). However, he took 
pressures in bars whereas they are in atmospheres in the 
original papers by Abraham et al. To check the effect of 




Pressure (bar) 

FIG. 9. Cube of the sound velocity as a function of pressure 
at 150 mK in helium 3 (circles) and helium 4 (squares). Linear 
fits (solid lines) give the spinodal pressure: Pa ~ —3. 14 bar in 
helium 3 and Ps ~ —9.65 bar in helium 4. 

this mistake on Pg , we performed the extrapolation with 
Eq. ( |Al| ) after the appropriate unit conversion. We used 
the same set of data points as did Maris, namely, the ones 
ranging from to 10 atm in helium 3 and from to 6 atm 
in helium 4. The linear fit of shown in Fig. ^ gives 
—3.1371 bar in helium 3 and —9.6456 bar in helium 4 for 
the spinodal pressure Ps- 

Since the experimental error bar is on c rather than on 
c^, we find it better to fit the data of Abraham et aZ.with 
the following formula: 

c^[b{P-Ps)]'/' . (A2) 

We now obtain Ps = —3.1534 bar in helium 3 and 
-9.6435 bar in helium 4, and b = 19.262 m'* s'^ kg^^ in 
helium 3 and 14.030 m^ s~^ kg~^ in helium 4. Note that 
we present the values with five digits in order to show 
the correction, although we think that the extrapolation 
does not attain this level of accuracy. 
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